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Abstract 

We study the "Casimir" energy of a minimahy coupled, real, massless scalar field outside a 
spherically symmetric background potential. We obtain a general expression for the null energy 
condition in d dimensions and explicit expressions for a perfectly reflecting spherical boundary in 
3+1 and 2+1 dimensions. In these cases, the null energy condition is always violated for radial 
motion and obeyed for azimuthal motion. Nevertheless, the averaged null energy condition is 
always obeyed. 
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I. INTRODUCTION 



There are several energy conditions in General Relativity which try to make a precise 
statement of the common sense notion that matter has positive energy density. The weak 
energy condition (WEC) states that no observer sees negative energy density, in other words 
that 



where T^i^ is the stress-energy tensor and is any timelike vector. If instead we consider 
null in Eq. (0), we get the null energy condition (NEC). 

The energy density between parallel conducting plates in the traditional Casimir effect ^ 
is negative, and thus provides one example of a system for which the WEC is violated. The 
pressure in the Casimir system is negative in the direction between the plates, so with 
pointing from one plate toward the other, the NEC is also violated. However, null vectors 
parallel to the plates have T^^V^V^ = because the positive pressure in that direction 
cancels the negative energy density. 

The averaged weak energy condition (AWEC) permits the violation of the WEC locally, 
but insists that when one integrates the WEC over a complete geodesic with tangent vector 
V^, the result will be non-negative. 



The averaged null energy condition (ANEC) is just Eq. (0) for null geodesies. 

One can imagine an observer at rest between the plates in the Casimir system, so that 
points only in the time direction. For that observer, the AWEC is violated. However, the 
Casimir system does not violate the ANEC, because geodesies parallel to the plates obey 
the NEC, while those that are not parallel eventually intersect the plates and pick up a large 
positive contribution from the plate material. 

The ANEC is in a certain sense the weakest of these conditions. If the ANEC is violated 
then the NEC must be violated somewhere, and so the WEC must also be violated, or else the 
NEC would hold by continuity. For a localized system in flat space, ANEC violation implies 
AWEC violation by the same continuity argument, but otherwise they are independent 
Nevertheless, the ANEC is sufficient to rule out many exotic phenomena, such as 
traversable wormholes j^, superluminal travel^ 0], and closed timelike curves and to 
prove singularity theorems 

It is very important to understand under what conditions the averaged null energy con- 
dition is violated, as Klinkhammer urged in . Klinkhammer found that for a free scalar 
field in Minkowski spacetime, the ANEC is satisfied. However, he found ANEC violations 
for the free scalar field in a fiat cylindrical topology, which corresponds to closing Minkowski 
space in one spatial direction. 

We study here a Quantum Field Theory (QFT) problem of a minimally coupled real 
massless scalar field outside a spherically symmetric background potential (i.e., a field with 
a position-dependent mass), in the calculational framework developed in [8|. In Sec. UTl we 
obtain a general expression for the null projection of the stress-energy tensor in arbitrary 



In this case, the relevant integral is not over a complete geodesic, but over the path to be traveled 
superluminally. 



T.^Vf'V > 



(1) 




(2) 
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dimension, followed by the results for 3+1 dimensions. In Sec. Illll we specialize to the case 
of a perfectly reflecting sphere in 3+1 dimensions. We find that the NEC is violated for 
radial motion and obeyed for azimuthal motion, whilst its validity for intermediate directions 
depends on the distance from the sphere. We then calculate the ANEC by integrating over 
a geodesic that passes outside the sphere. Although the NEC is violated far from the sphere 
where the motion is primarily radial, the points closer to the sphere dominate, and as a 
result the ANEC is always obeyed. 

In Sec. II VI we briefly report our calculations of the NEC and ANEC outside a perfectly 
reflecting circle in 2+1 dimensions. The results are very similar to those for 3+1 dimensions, 
and once again we find that although the NEC may be violated, the ANEC is satisfied. 

II. NULL ENERGY CONDITION FOR A SCALAR FIELD 
A. General case 

In this section we find a general expression for the NEC, for the case of a minimally 
coupled, real scalar field outside a potential. The Lagrangian is 

L = \ {d,<^d^<^ - (3) 

and the stress-energy tensor is 

T^. = d^<\>dA + \]^v ^{x) - aV^A^] . (4) 
For a null vector, ri^yV^V^ = 0, so we have 

T.^V^V = . (5) 

For a static system^, Toj = 0. If we further choose coordinates in which Tij is diagonal, and 
V = (l,v), then 

T^uV^V = ^' + Yl i^A<pf . (6) 



B. Spherical symmetry 

We are going to calculate the NEC for the scalar field outside a potential that is spherically 
symmetric in m spatial dimensions. Using spherical coordinates and choosing V to lie in 
the equatorial plane, we find 

T^.V^'V = 02 + {vrdAf + [vj^d^'i^ (7) 

where Vr is the component of the velocity in the radial direction and is the azimuthal 
component. Since we are taking = 1, v'i + = I- 



^ Greek letters run over all indices and Latin ones run over spacelike indices only. 
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Decomposing the quantum field in terms of modes gives 



0(r,^],t) = (8) 

Jo V2.7[UJ 

wliere u = k since tlie field is massless, and the sum over i gives the partial wave expansion 
in the m spatial dimensions. 

The wavefunctions ipl{r) are the eigenstates of the time- independent radial Schrodinger 
equation 

(f m - 1 d i(i + m-2) ^^ x , 2 -/'z x 

- ^— - + ^ + Vir)j ^iir) = e^iir) . (9) 

In general, the solutions to Eq. ^ comprise both bound and scattering states, but here we 
will consider only repulsive potentials, so there are no bound states. 

The wavefunctions and creation and annihilation operators are normalized as follows: 



27r^ 



= (lOa) 

r""~^^i{ry^t{r) dr = Tx5{k - k') (10b) 

[afn^'^l = [<%4"1=0 (10c) 

= ^{k - k')5ee'5e,e', . (lOd) 

Using these expressions, the vacuum expectation of the time derivative term in Eq. ((Tj) 
is found to be 

^E^r/°°fc.(|^f(r)p-|^f(°)(r)p) . (11) 



where we have used dimensional regularization |8|, and have subtracted the free wavefunc- 
tions ipl^'^\r). No other counterterms are necessary, because we are considering only lo- 
cations outside the potential. The factor is the degeneracy in each partial wave 



T{m + i-2) 

"r(m-i)r(£+i)^'" + ^^ ^> ^^^^ 



with D^ = l. 



The vacuum expectation of the radial term of Eq. ((7j) gives 

k^E^r (\dri^i(^)\' - \dri^f\^)\') ■ (13) 

In the azimuthal term of Eq. ((7j) we have to differentiate the spherical harmonic, Y^™ (f2). In 
2+1 dimensions, Y^™(f2) oc e"^^""^. The same result holds for 3+1 dimensions when we consider 
the specific case of geodesies in the equatorial plane. Thus, the vacuum expectation of the 
azimuthal term gives 

kE^r r-—,mir)\'-\^i^'\r)\^) (14) 
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where 



2£2 m = 2 

1)(2^+ l)/2 m = 3. 



(15) 



Adding Eqs. (HH), (0 and dHj), we find 



1 ^ dk 



UJ 



Outside any spherically symmetric potential, the wave functions are given by 

(r) = yAUO^ [e''''Hi'\kr) + Hi'\kr)] 



(16) 



(17) 



where 6e is the scattering phase shift in the quantum mechanical problem with the same 
potential, z/ = m/2 — and the normalization factor 



r(m/2) 



The free wavefunctions 'ipl^^\r) are given by Eq. p7|l . with 6e = 0. Thus we find 

Mr)\' - |^r(r)P = NUr)k [{e'^'^ - l) Hi'\krr + {e~'^'' - l) Hl'^krY] 
The radial derivative in Eq. (fT^ is 



where 



and likewise for Hy . Therefore 

\dri,i{r)\' -\d,.i,T\r)\' = N^{r)e \{i 



2i5e 



1) Hl'\krf + {e-^''^ - 1) Hj^\kr 



(18) 

(19) 

(20) 
(21) 

(22) 



We show in the appendix that the second term of Eq. (jl9|) is just the first term with the 
replacement k — > —k + ie and likewise for Eq. (j22I)- Thus we can drop the second terms and 
extend the range of integration over k to — oo, with the understanding that k is to be taken 
above any branch cut on the negative real axis. Then 



dk 



X 



{e^'''-l)NUr)k 



Hi'\krf + Dl 



Hj^\kr 



(23) 



Following the methods used in we now convert this expression to a contour integral 
which we close in the upper half plane. In general, 6i{k) will not be well behaved off the real 
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axis, but in the present case of a potential with compact support there will be no difficulty. 
The only contribution to the integral comes from the branch cut along the imaginary k axis. 
To the right 00 = Vk^ = k, but to the left uj = —k, so with k = in, we obtain 



TT 



POO 

dK I {^''^^'^^ - 1) 

JO 



(24) 



DT-a 



Using Eq. da} and H^y\ix) = {2/n)i-^''+^^K^{x), we get 



4iVn^(r; 



[Krf (25) 



Kr 



for the general case of a minimally coupled, massless, scalar field outside a spherically 
symmetric potential. 

For m = 3, we find u = / + 1/2, A/'„(r) = l/(16r), = 2£ + 1, = £(£+ 1)(2£+ l)/3, 
and 



1 r°° 

T^.V^V = —y- V(-l)' / (e2*^^(-) - 1) 



(26) 



2£ + 1 - 



' Kf+i/2(Kr) 



2K^r2 



KT 

In Sec. IIVI we will summarize the results for 2+1 dimensions, which closely mimic our 
findings for 3+1 dimensions. 

III. PERFECTLY REFLECTING BOUNDARY CONDITIONS 
A. Phase shifts 



For a hard sphere of radius a, the scattering phase shift is 

H^u\ka) 



Hi^\ka) 

so that the wave function vanishes at r = a. Thus 

2J^{ka) 



Hi'\ka) 



(27) 



(28) 
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FIG. 1: We consider the null energy condition at distance r from the center of a sphere for motion 
in the direction V which makes an angle a with the radial direction. 



and 



2v 



IT 



Inserting Eq. ()29|) into Eq. ()25p. we find for general m, 



7r2 ^ ./o 



and for m = 3, 



Ka) 



1 - 



(29) 



(30) 



Kr 



(31) 



B. NEC numeric results 3+1 dimensions 



Equation (j!?T|l gives T^jjV^V^ in terms of the radial distance from the sphere and the 
components of velocity in the radial and azimuthal directions. We can rewrite it in terms of 
the angle a that any arbitrarily chosen velocity vector makes with the radial direction, as 
depicted in Fig. [TJ so that t>^ = sin^ a, and = cos^ a. We find 



£{i + 1) sin^ a 



(32) 



cos a I — KJkt) — Kyj^iinr) 
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NEC VIOLATION 



1.5 
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FIG. 2: Region of NEC violation (shown shaded) in 3+1 dimensions. The parameter r is the radius 
in units of a, and a is the angle between V and the radial direction. 
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FIG. 3: We integrate the NEC along a line parallel to the x-axis with impact parameter h. 

We have used Mathematica to determine numerically where the NEC is and is not obeyed. 
Figure 121 shows the values for which the NEC is violated, as a function of radial distance r 
and angle a. Radial vectors always violate the NEC, while azimuthal vectors always obey 
it. The dividing line depends on the distance from the sphere. 

C. ANEC in 3+1 dimensions 



We obtain the ANEC by integrating the NEC along a complete null geodesic. We consider 
an observer moving along a geodesic parallel to the x-axis and passing by the spherical 
boundary with an impact parameter 6, as depicted in figure El We find 



/CXJ -1 /'OO /»00 / T I 



IJko) 



na] 



2I 
K - 



(33) 



X 



1 - 



KJnr) 



Kr 



K^{Kr) - K^+i(Kr) 
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FIG. 4: Contributions to the ANEC integral from the first three partial waves = solid, I = 1 
dashed and I = 2 dotted) in 3+1 dimensions vs impact parameter 6, with a = 1. 
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FIG. 5: ANEC integral in 3+1 dimensions vs. impact parameter h, with a = 1. 

where = + 6^, t;^ = sin^ a = + 6^), and = cos^ a = /{x^ + 6^). 

We have found numerically that the ANEC is indeed obeyed. The sum over partial waves 
converges rapidly. The £ = partial wave always obeys the NEC, but the other partial waves 
do have negative contributions for certain impact parameters b, as shown in Fig. 0] However 
these negative contributions are always more than compensated for by positive contributions 
from lower as can be seen in Fig. El which shows the contribution including all partial 
waves. 



IV. NEC AND ANEC IN 2+1 DIMENSIONS 

The results obtained for 2+1 dimensions closely follow those for 3+1 dimensions. For 
m = 2, we find 1/ = / and Nm{r) = 1/8. With / = 0, we use = 1, = and Eq. 
becomes 

T.^V^V = — dn {e^^'^('^^ - 1) [Mnrf - v^K^inrf] . (34) 
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NEC VIOLATION 
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FIG. 6: Region of NEC violation (shown shaded) in 2+1 dimensions. 



For / > 0, we use Dj = 2, C| = 2/^ and we get 



7^0 



^0 



With the perfectly reflecting circular boundary condition for / = we obtain 

1 r . loi^a) 



27r2 



dn 

KQ{Ka 



and for / > we find 



-y 

7r2 ^ 



Kpina] 



KT 



(35) 



(36) 



(37) 



Figure ini shows the values for which the 2+1 dimensional NEC is violated, as a function 
of radial distance r and angle a. Once again, radial motion violates the NEC. 
As in 3+1 dimensions, the ANEC is always obeyed, as shown in Fig. [7| 

V. CONCLUSIONS 



We have studied the problem of a minimally coupled, real, massless scalar field outside a 
spherically symmetric background potential, within the calculational framework developed 
in j^. We obtained a general expression for the null energy condition in d dimensions. We 
calculated the NEC for the specific case of a perfectly refiecting sphere and showed that the 
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FIG. 7: ANEC integral in 2+1 dimensions vs. impact parameter b, with a = 1. 

NEC is obeyed at some points in space and violated at others, depending on the velocity of 
the observer. 

We calculated the averaged null energy condition for 3+1 and 2+1 dimensions by integrat- 
ing the NEC over a complete null geodesic. Although there is a range of impact parameters 
for which partial waves with i > contribute negatively to the ANEC, the ANEC is always 
obeyed when summing over all the partial waves. 

Although we have done explicit calculations only for perfectly reflecting boundary con- 
ditions, we conjecture that the ANEC is satisfied for all geodesies which pass outside any 
localized background potential. 
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APPENDIX A: ANALYTIC PROPERTIES OF THE PHASE SHIFT 

In this appendix we show that, for positive, real k, 

(_A;) [e^^^d-k) _ 1) Hl'\-krf = k (e'^^^^W - l) Hl^\krf (Al) 

where —k is taken in the upper half plane, by using the techniques of 0. From Eqs. (22) 
and (21) of 0, and using the notation of that paper, we find the Green's function, 

-TTi{iy-l/2) 

Gi{r, r, k) = Mk, r)fi{k, r) . (A2) 

k 

Outside the potential, the Jost solution is just the free Jost solution, 

f,{k,r) = e^^-'^^H^'\kr) (A3) 
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while the physical scattering solution can be written in terms of the scattering phase shift, 

Mk, r) = ^ [e'^'^^'^Hi'\kr) + Hi'\kr)] . (A4) 



Thus 



Gi{r,r,k) = '^Hl^\kr)[e^'^'^^^Hl^\kr) + Hi^\kr)] (A5) 
= ^ [{e'^''^"^ - 1) Hi'\krf + 2Mkr)Hi'\kr)] . 

Now G has the property G{r, r, —k) = G{r, r, k)*, and thus we find 

(e^'^^(-''') - 1) Hl^)(-kr)^+2J,{-kr)Hl^\-kr) = - (e~^'^'^^^ - l) Hl^\kry-2J,{kr)Hl^\kr) . 

(A6) 

The second terms on the two sides are equal, so the first terms must be equal also, which 
proves Eq. ()A1|) . 

As A; — i> —k, H^\krY transforms in the same way as so 

(_A;)3 [e^^'d-^) - 1) H^;)(-krf = P [e-^^^^^^ - l) H'i^\krf (A7) 

as well. 
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